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ABSTRACT 


A  comparison  is  made  between  theoretical  and  experimental  results 
for  cooperative  scattering  between  two  spheres.  The  overall  agreement 
between  theory  and  experiment  is  quite  good.  Also  a  large  side  scatter¬ 
ing  resonance  which  was  measured  to  be  44  times  larger  than  that  due  to 
a  single  sphere  was  calculated  to  be  actually  47.6  times  larger. 


I 


In  a  recent  publication,  Wang  et  al.l  presented  some  extremely 
interesting  measurements  from  the  cooperative  scattering  from  two 
spheres  using  microwave  analog  techniques.  The  results  they  presented 
in  the  above  reference  were  only  a  small  part  of  a  more  comprehensive 
study2.  One  of  the  more  interesting  discoveries  in  this  study  was  the 
enormous  enhancement  of  the  scattered  intensity  component  perpendicular 
to  the  scattering  plane,  at  right  angles  to  the  incident  beam  for  two 
identical  spheres  in  contact  when  the  axis  of  symmetry  of  the  two 
spheres  bisected  the  scattering  angle. 

It  is  the  purpose  of  this  short  paper  to  demonstrate  that  their  re¬ 
sults  are  borne  out  beautifully  by  theoretical  calculations.  We  have 
used  the  powerful  technique  presented  by  Bruning  and  lo3>4  to  con¬ 
struct  a  very  complex  computer  program  to  perform  the  calculations.  It 
should  be  noted  that  several  errors  (probably  typographical)  were  disco¬ 
vered  in  their  equations  and  virtually  all  of  them  had  to  be  rederived. 
In  a  future  publication  we  will  present  a  corrected  set  of  equations 
along  with  several  improvements  for  other  researchers  to  use.  The  pro¬ 
gram  was  written  in  Fortran  77  employing  complex,  double  precision 
arithmetic.  The  program  was  subjected  to  many  rigorous  tests.  In  Fig. 

1  we  show  the  results  for  extinction  coefficient  calculations  as  a  func¬ 
tion  of  distance  of  separation  of  centers  (kd=2wd/A).  The  case  shown  is 
for  a  size  parameter  x=0.9283  and  refractive  index  N=1.54.  If  this  case 
Is  compared  to  a  similar  case,  calculated  by  Kattawar  and  Humphreys^ 
(Fig.  6  of  that  reference)  using  the  point  dipole  approximation,  which 
Is  a  completely  Independent  method  and  program,  it  will  be  seen  that  the 
agreement  is  remarkably  good. 

We  next  calculated  the  following  quantities  to  compare  with  the 

2 

measurements  of  Schuerman  and  Wang  : 


1)  |S(0)|  -  modulus  of  complex  scattering  amplitude  In  the  forward 

direction 

2)  P  I"i[S(0)] 

k^G 

3)  Q  -  X-  Re[S(0)3 

k^G 

4)  <fr(0)  -  phase  angle  defined  as  the  angle  between  S(0)  and  the 

P  axis 

where  G=2»a2  where  a  is  the  sphere  radius  and  k  is  the  wavenumber 
2t r/X.  Table  I  gives  the  comparison  between  theory  and  experiment  for 
three  different  size  parameters  and  refractive  indices  when  the  spheres 
are  In  contact.  Also  for  ease  of  visualization  we  have  also  presented 
the  P  versus  Q  plots  in  Figs.  2a-2c.  There  are  several  noteworthy  fea¬ 
tures  in  the  comparisons.  First,  the  agreement  between  theory  and  ex¬ 
periment  on  the  phase  angles  is  within  a  few  percent  in  all  but  one  mea¬ 
surement  (x-3. 752,a=0o)  where  it  reaches  8%.  Secondly,  the  worst  disa¬ 
greement  in  all  parameters  usually  occurs  at  a=0°,  i.e.,  for  end  on  in¬ 
cidence  which  is  the  region  most  difficult  to  measure.  pinally,  in  the 
P-Q  plots  of  Fig.  2,  the  theoretical  results  always  converge  closer  to 
the  non-interacting  sphere  limit  (labelled  NIS  on  the  figures)  for  a=90° 
(broadside  incidence)  than  the  measurements.  This  fact  is  borne  out 
more  realistically  in  Fig.  2d  where  the  spheres  are  separated  by  3.71 
diameters. 

The  truly  exciting  results  are  presented  in  Figs.  3a  and  3b.  In 
Fig.  3a  we  show  the  intensity  component  perpendicular  to  the  scattering 
plane  for  a  90°  scattering  angle  (Ir(90°))  as  a  function  of  particle 
orientation  angle  a.  If  this  case  is  compared  with  the  experimental  re- 


suits  of  Schuerman  and  Wang2  (their  Fig.  12A)  the  only  major  differ¬ 
ence  is  the  maximum  value  where  they  measure  37.5  and  the  calculations 
give  43.6.  Also  shown  on  this  figure  is  the  sensitivity  of  this  peak  to 
particle  separation.  When  the  spheres  are  separated  by  1.5  sphere  dia¬ 
meters  (kd-14.034)  the  maximum  almost  becomes  a  minimum,  decreasing  more 
than  a  factor  of  ten.  We  also  calculated  the  parallel  component  of  the 
radiance  when  the  electric  vector  is  parallel  to  the  scattering  plane. 

It  also  shows  a  maximum  at  o=45°  but  it  is  roughly  a  factor  of  two  smal¬ 
ler.  To  see  if  another  maximum  occurred  in  the  orthogonal  position  we 
rotated  the  spheres  in  the  opposite  direction  so  that  at  a=45°  the  axis 
of  symmetry  of  the  spheres  was  perpendicular  to  the  bisectrix.  Another 
maximum  occured  but  it  was  roughly  five  times  smaller  than  the  maximum 
which  occurred  when  the  axis  of  symmetry  was  in  the  specular  position. 

In  Fig.  3b  we  show  a  similar  plot  except  the  size  is  changed  and  the  re¬ 
fractive  index  is  now  complex. This  case  can  be  compared  with  Fig.  2  of 

Wang  et  all.  Again  the  overall  agreement  is  extremely  good  but  as  in 

• 

the  previous  case  we  calculate  a  higher  maximum,  namely  50  compared  to 
46.2  or  a  factor  of  47.6  times  larger  than  the  radiance  for  a  single 
sphere.  If  the  spheres  are  moved  slightly  apart  (1.04  sphere  diameters) 
the  maximum  falls  to  44.5.  If  they  are  two  diameters  apart  the  maximum 
becomes  a  minimum  but  smaller  secondary  maxima  start  to  appear.  These 
results  show  that  the  combination  of  multiple  scattering  of  the  near 
fields  of  the  two  spheres  with  the  direct  scattered  fields  have  just  the 
right  phase  behavior  to  produce  the  enhancement. 

The  overall  agreement  between  theory  and  experiment  for  two  spheres 
cooperative  scattering  is  quite  good.  In  a  future  publication  we  will 
explore  the  enhancement  phenomena  in  greater  detail  as  well  as  study  the 
complete  Mueller  matrix. 


This  research  was  partially  supported  by  National  Science  Founda¬ 
tion  grant  No.  ATM-8204191  and  the  Office  of  Naval  Research  through  con¬ 
tract  N00014-80-C-0113. 
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FIGURE  CAPTIONS 


Fig.  1.  Extinction  coefficient  calculations  for  two  spheres  as  a 

function  of  distance  of  separation  (kd).  The  spheres  have  a 
size  parameter  x=0.9283  and  a  refractive  index  N=1 . 54.  The 
case  o=0°  corresponds  to  end-on  incidence  while  a=90°  corres¬ 
ponds  to  broadside  incidence.  The  subscript  I  refers  to  the 
incident  beam  polarized  in  the  scattering  plane  while  r  is 
for  incident  polarization  perpendicular  to  the  scattering 
plane. 

Fig.  2.  (a-c)  P-Q  plots  for  the  comparison  between  theory  and  experi¬ 

ment  for  different  sizes  and  refractive  indices  for  two  spheres 
in  contact  and  for  orientations  from  end-on  incidence  (0°)  to 
broadside  indicence  (90°).  The  point  labelled  NIS  is  the  P-Q 
value  obtained  for  non-interacting  spheres.  Curve  d  is  for  a 
case  of  larger  separation. 

Fig.  3.  (a-b)  Scattered  radiance  at  90°  perpendicular  to  the  scatter¬ 

ing  plane  (Ir  (90°))  as  a  function  of  orientation  angle  a. 

Note  at  o=45°  the  symmetry  axis  of  the  particle  bisects  the 
scattering  angle  corresponding  to  specular  reflection. 
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Figure  3 


